METHOD FOR THE SYNCHRONIZATION OF DATA, SUCH AS 
DISTRIBUTED DATA, TAKING ACCOUNT OF CLOCK DRIFTS AND 
INACCURACIES. 

The invention relates to a method for determining a 
5 correspondence between the time indicated by the internal clocks of a 
plurality of machines connected to a network, such as sensors or computere, 
and also to a method for synchronizing data originating from said machines. 
These two types of process are intimately linked, the second constituting a 
direct application of the first. 
^° 'n this context, 'synchronizing data' signifies ordering them 

chronologically as a function of their time of acquisition; In a more restricted 
sense, it also signifies 'date stamping' these same data values with respect to 
a unique temporal reference. On the other hand, 'synchronizing clocks' 
signifies forcing the latter to indicate the same time at a given moment. 

A method of synchronizing data is necessary, for example, in 
the field of the validation of automobile driving aid devices. In order to study 
the effect of these devices on the behavior of a driver, test automobiles are 
equipped with various sensors connected In a network and allowing 
information on the vehicle (average speed, acceleration...), on the interaction 
2 0 with other vehicles (distance and relative speed. . . ) and on the behavior of the 
driver (reaction time, direction of observation...) to be simultaneously 
acquired. The data values acquired by these various sensors and those 
recorded by the driving aid devices must be synchronized in order to be 
exploitable by an ergonomist, and the precision of the synchronization must 
25 be better than one millisecond, and preferably better than 100 ^is, for data 
acquisition rates able to reach 1 kHz. At the same time, it Is desirable to use 
standard equipment devices and for the system to be very flexible, allowing 
the acquisition network to be reconfigured during use, and to be very robust 
against failures and connection problems. 
30 The most obvious solution to the problem of the 

synchronization of the data originating from a plurality of devices consists in 
synchronizing all the clocks. Each device assigns a 'date stamp' to the data 
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values it records and, since the clocl<s are synchronized between one 
another, all the date stamps are directly comparable. In this context, the 
problem of determining a correspondence between the local times is not 
posed, because all these times substantially coincide. 
5 The synchronization of the clocks may be performed in a 

hardware or software manner: a general introduction to these techniques is 
provided in the article by J.A. Fonseca and P. Fonseca, 'An overview of clock 
synchronization solutions', 4*'' IFAC International Symposium, Sicica, 
Argentina, 2000. 

10 The hardware solutions are based on the use of an external 

clock signal, common to all the devices. In this way, a very high precision can 
be obtained, for example by the use of an atomic clock, such as that of the 
GPS system. On the other hand, these solutions have the drawback of 
requiring dedicated equipment: devices that are not designed to use an 

15 external clock signal cannot be connected to the network. 

The software solutions comprise the use of software clocks 
(variables stored in memory) controlled by the 'hardware' clock of each device 
connected to the network. These software clocks are synchronized in such a 
manner as to indicate approximately the same time at any moment. Among 

2 0 these solutions, two classes can be further differentiated: 

- either a single synchronization is performed at the 
beginning of the session using the network, and it is assumed that the drifts of 
the clocks can be neglected; 

- or, re-synchronizations are also performed during use; 
25 this is the technique adopted, for example, by the protocol NTP (Network 

Time Protocol). 

As regards the first class of solutions, the assumption of 
negligible drifts is difficult to verify, at least for acquisition times of several 
hours and when commercially available computers are used: the clocks of the 
30 latter exhibit drifts that may reach 300 |js/s, or over one second per hour. 
Replacing these clocks by components of better quality equates to a 
hardware solution that is costly and not very flexible. 
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The second class of solutions has a difficulty associated with 
the monotonicity of the clocks. If, for example, a clock has been gaining 
relative to the reference time, it will be turned back during the re- 
synchronization. If a first data value has been recorded immediately prior to 
5 the moment of re-synchronization and a second data value immediately 
aftenwards, the time stamp of the second recording pre-dates that of the first 
and the temporal order of the data is reversed. Moreover, the known 
algorithms that impose the conservation of the monotonicity of time exhibit a 
relatively slow convergence. 

10 In any case, basing the synchronization of the data on the 

synchronization of the clocks, or on the use of a single clock, poses limits to 
the possibilities for reconfiguring the network of devices. Here is a concrete 
example of a case where the limits of the known techniques of the prior art 
can be seen. Considering two machines, X and Y, connected to the same 

15 network and whose clocks are kept synchronous by a hardware or software 
method. At a given moment, the machine Y is disconnected from the network, 
purposely or because of a connection problem, but continues to record data; 
finally, the machine Y is reconnected after a certain lapse of time. It can easily 
be understood that it is not possible to order the data acquired during the 

20 period of disconnection chronologically, because the synchronization between 
the clocks of the machines X and Y is irretrievably lost. For this reason, these 
methods of synchronization are very sensitive to potential network connection 
problems. In addition, if during the period of disconnection the machine Y had 
been connected to a second network also comprising a machine Z, there 

25 would be no means of synchronizing the data recorded by X with those 
recorded by Z. The reconnection of the machine Y is also problematic, 
especially if it is desired that the monotonicity of its clock be conserved. For 
this reason, these known synchronization methods pose limits to the 
possibilities for 'hot' reconfiguration (during operation) of the networks to 

3 0 which they are applied. 

Another technique known from the prior art is to use temporal 
markers in order to establish an approximate con-espondence between the 



device local times whose internal clocks are not synchronized. In this case, 
the synchronization of the data is most often done retro-actively, in other 
words after the end of the recording session, and for this reason it is 
sometimes referred to as 're-synchronization' of the data. 
5 Such a case frequently arises in the field of multimedia, where 

it is for example necessary to synchronize images, acquired by a digital 
camera, with sounds, acquired by an external microphone. In the field of 
multimedia, the demands in temis of precision are however fairly modest, 
errors up to around 33 milliseconds (ms) being undetectable by human 

10 beings. Several data synchronization techniques used in multimedia 
applications are described in the article by G. Blakowski and R. Steinmetz, 'A 
media synchronization survey: reference model, specification, and case 
studies', IEEE J. Selected Areas Commun. 141 (1996), pages 5 - 35. 

These techniques are robust and allow a wide flexibility, but 

15 their precision is very limited and insufficient for many applications. Moreover, 
it is not possible to determine with certainty an upper limit for the amplitude of 
the synchronization errors committed. 

A subject of the present invention is a method for establishing 
a correspondence between the local times of two or more machines whose 

2 0 clocks are not synchronized. 

Another subject of the present invention is such a method, 
with an improved precision with respect to the prior art. 

A further subject of the present invention is such a method, 
with a known precision modeled by an interval. 

25 A further subject of the present invention is such a method, 

which is robust against connection problems of the network to which the 
machines, whose local times are to be made to correspond, are connected. 

A further subject of the present invention is such a method, 
which allows a dynamic reconfiguration of such a network. 

30 A further subject of the present invention is such a method, 

which allows only standard equipment to be used. 



A further subject of the present invention is a method for 
synchronizing the data recorded by two or more machines, or generated by 
two or more devices, which does not require synchronization of the clocks of 
said devices or machines, and which is based on the establishment of a 
5 correspondence between the local times of two or more machines whose 
clocks are not synchronized. 

Other objects of the present invention consist in providing 
such a method with an improved precision relative to the prior art, with a 
known precision, that is robust against problems of connections of the 
0 network to which the machines, whose local times are to be made to 
correspond, are connected, and/or which allows only standard equipment to 
be used. 

At least one of the aforementioned objects is achieved by 
means of a method for establishing a correspondence by intervals between 

5 the time indicated by a first monotonic clock and the time indicated by a 
second clock, also monotonic, characterized in that it operates, over at least 
one temporal range, a first temporal reference common to the first and to the 
second clock and monotonic over said range, and in that said method 
comprises the steps of: 

3 a) determining a first temporal interval bounded by a 

first pair of time values of the first clock and belonging to a first temporal 
range over which said first temporal reference common to the first and to the 
second clock exists; 

b) determining a second temporal interval bounded by 
5 a second pair of time values of the first clock and belonging to a second 

temporal range over which said first temporal reference common to the first 
and to the second clock exists; 

c) determining, using the common temporal 
reference, a third temporal interval, bounded by a first pair of time values of 
the second clock, and containing the first temporal interval; 
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d) determining, by means of the common temporal 
reference, a fourth temporal interval, bounded by a second pair of time values 
of the second clock, and containing said second temporal interval; 

e) for any given fifth temporal interval bounded by a 
5 third pair of time values of the first clock, calculating a sixth temporal interval, 

bounded by a third pair of time values of the second clock, and containing 
said fifth temporal Interval, the calculation being performed by interpolation or 
extrapolation using said first, second, third and fourth intervals. 

In one particular embodiment of the invention, the step c) 
10 comprises the operations of: 

c1) determining a seventh and of an eighth temporal 
interval, bounded by a fourth and fifth pair of time values of the second clock, 
respectively, and belonging to the first temporal range, such that the first 
temporal interval falls in the range between said seventh and eighth temporal 
15 intervals; 

c2) determining a first, second and third value of the first 
common temporal reference, included in said first, seventh and eighth 
temporal intervals, respectively; 

c3) calculating, by interpolation, said third temporal interval, 
20 using said first, seventh and eighth temporal intervals and said first, second 
and third values of the first common temporal reference; 

and the step d) comprises the operations of: 

d1) determining a ninth and of a tenth temporal interval, 
bounded by a sixth and seventh pair of time values of the second clock, 
25 respectively, and belonging to the second temporal range, such that said 
second temporal interval falls in the range between said ninth and tenth 
intervals; 

d2) determining a fourth, fifth and sixth value of the first 
common temporal reference, included in said second, ninth and tenth 
30 temporal intervals, respectively; 



d3) calculating, by interpolation, said fourth temporal 
intervals, using said second, ninth and tenth temporal Intervals and said 
fourth, fifth and sixth values of the first common temporal reference. 

Preferably, the operation c3) is carried out by replacing said 
5 first, second and third values of the common temporal reference by temporal 
intervals whose width is equal to the discretization of the firet common 
temporal reference over the first temporal range, and the operation d3) is 
carried out by replacing said fourth, fifth and sixth values of the common 
temporal reference by temporal intervals whose width is equal to the 
10 discretization of the first common temporal reference over the second 
temporal range. 

Advantageously, during the temporal range or ranges over 
which a first common temporal reference exists, a first reading of the first 
clock is recorded several times, followed by a reading of the first common 

15 temporal reference and, subsequently, by a second reading of the first clock, 
and independently, a first reading of the second clock is recorded, also 
several times, followed by a reading of the first common temporal reference 
and then by a second reading of the second clock, and the operations c1), 
c2), c3), d1), d2) and d3) are performed using these recordings. 

20 In another particular embodiment of the invention, the step c) 

comprises the operations of: 

c1) determining a seventh and an eighth time value of the 
second clock belonging to the first temporal range, such that said first 
temporal interval falls in the range between said seventh and eighth values; 

2 5 c2) determining a first, second and third interval of values of 

the first common temporal reference, comprising said first temporal interval 
and said seventh and eighth time values of the second clock, respectively; 

c3) calculating, by interpolation, said third temporal interval, 
using said first interval of time values of the first clock, said seventh and 

30 eighth time values of the second clock and said first, second and third 
intervals of values of the first common temporal reference; 

and the step d) comprises the operations of: 
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d1) determining nintfi and tenth time values of the second 
clocl< belonging to the second temporal range, such that said second temporal 
interval falls in the range between said ninth and tenth values; 

d2) determining fourth, fifth and sixth intervals of values of 
5 the first common temporal reference, comprising said second temporal 
interval and said ninth and tenth time values of the second clock; 

d3) calculating, by interpolation, said fourth temporal 
interval, using said second interval of time values of the first clock, said ninth 
and tenth time values of the second clock and said fourth, fifth and sixth 
10 intervals of values of the first common temporal reference. 

Preferably, the operations c3) and d3) are carried out by 
replacing said seventh, eighth, ninth and tenth time values of the second clock 
by temporal intervals whose width is equal to the discretization of the time of 
the second clock. 

15 Advantageously, during the temporal range or ranges over 

which a first common temporal reference exists, a first reading of the first 
common temporal reference is recorded several times, followed by a reading 
of the first clock and, subsequently, by a second reading of the first common 
temporal reference, and independently, a first reading of the first common 

20 temporal reference is recorded, also several times, followed by a reading of 
the second clock and then by a second reading of the first common temporal 
reference, and the operations c1), c2), c3), d1), d2) and d3) are performed 
using these recordings. 

In particular embodiments of the invention, the 

25 aforementioned calculations by interpolation or extrapolation are, more 
precisely, interpolations or extrapolations that are linear or linear by intervals. 

Advantageously, when the first temporal reference common to 
the first and to the second clock exists over at least two separate temporal 
ranges and can comprise a rupture in monotonicity from one temporal range 

30 to the other, a second monotonic common temporal reference is used in order 
to resolve the ambiguities resulting from the non-monotonicity of the first 
common temporal reference. 
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The invention also relates to a method for synchronizing data 
recorded and date stamped by a first machine, having a first clock, with 
respect to the local time of a second machine, having a second clock, 
characterized in that the date stamping is carried out by associating with each 
5 data value recorded by the first machine a fifth temporal Interval, bounded by 
a third pair of time values of the first clock, and in that the synchronization is 
performed by determining, using a method such as that described 
hereinabove, a sixth temporal interval, bounded by a third pair of time values 
of the second clock and containing said fifth temporal interval. 
10 In a particular embodiment of such a data synchronization 

method, the first common temporal reference is supplied by the clock of a 
synchronous bus which connects, at least temporarily, said first and second 
machines. 

In a particular embodiment of the invention, such a data 
15 synchronization method may be broken down into a first sub-process for 
recording clock readings, such as that described hereinabove, carried out 
locally by each machine whose data it is desired to synchronize, and a second 
process for the synchronization itself, effected by a single machine toward 
which all the data has been transferred after the termination of the data 
2 0 recording session. 

Other features, details and advantages of the Invention will 
become apparent upon reading the description presented with reference to 
the appended drawings, given by way of example and in which: 

figures 1A and 1B illustrate a manner in which to proceed for 
25 the determination of a correspondence between the times indicated by two 
different and non-synchronized clocks; 

figures 2, 3A - 3L, 4 and 5, together with tables Tab.A, Tab.B 
and Tab.C, illustrate by a concrete example a method forming one 
embodiment of the invention; and 
30 figure 6 shows a flow diagram of this method. 

Before continuing with the description, it is appropriate to 
define precisely certain terms and notations used in the following description. 
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By 'true time' is understood physical time, whicli cannot be 
precisely known and of which the clocks only give an approximation; in this 
document, 't' indicates a true time value. 

By 'clock' is understood a digital clock, formed by an oscillator 
coupled to a counter. At each cycle of the oscillator, the counter is 
incremented by a discrete quantity 6 (discretization or quantization). 
Consequently, the plot of the time measured by a clock with respect to true 
time takes the form of a staircase, with steps of height 6. The counter of a 
clock necessarily has a finite capacity and is reset to zero when this is 
exceeded. Clocks can nevertheless have a capacity that is sufficiently large 
for this problem to be neglected. 

A 'synchronous bus' is a bus over which a clock signal is 

broadcast. 

By 'local time' of a machine is understood the time measured 
by the clock of this machine. A local time value is indicated by T. 

By 'date stamp' of a data value recorded by a machine is 
understood either the time read on the clock of the machine at the moment 
the data value is recorded, or the interval included between two readings of 
this clock, carried out before and after said recording. Any ambiguity will be 
removed by the context and by the notation: thus, the date stamp indicated by 
'[T]' is, in fact, the interval included between land f , with T<f . In reality, 
any date may be considered as an interval: even if it is known with certainty 
that it was the value indicated by the clock of the machine at the moment a 
data value was recorded, there remains an indetermination 6 due to the 
discretization of said clock. 

By 'monotonicity' is understood the property of a clock such 
that if ti<t2, then T(^)< T{t2). 

By 'offset' is understood the difference between the time 
indicated by two clocks at a given moment, or between the time on a clock 
and true time. The offset between two times defined by means of intervals is 
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also an interval: if [Ta] = [T^ , Ta] and [Tb] = [Ib ,Tb], then the offset allowing 
the passage from [Tb] to [Ta] is equal to: 

[0ffAB] = [0ff^, ^] = [Ia-Tb, Ta-Ib] (1) 
such that [Ta] = [Tb] + [offAs]. 
5 By 'drift' between two clocks, or one clock and true time, is 

understood the derivative of the offset with respect to local time of one of the 
clocks or of true time. Unless othenA/lse stated, in this document it is 
considered that, for each pair of clocks, the drift is constant over time and 
that, consequently, the offset is a linear function of time (linear drift 
10 hypothesis). In this case, the drift can be calculated by knowing two offsets 
determined at different dates. If the offsets and/or the date stamps are 
intervals, the drift is also an interval: 



[driftAB] = [ drift^B . drift AB ] = ^"-^ 1 



in the sense of the calculation by intervals. 

The exponents 1 and 2 refer to the first and to the second 
date at which an offset is determined. 

By 'synchronization by intervals of the data of the machine B 
20 with respect to the machine A' is understood the determination, for each data 
value recorded by the machine B with a date stamp [Te], of an interval [Tab] of 
the local time of the machine A such that precedes the recording of said 

data value and T^b follows it. For the sake of concision, [Tab] is said to be the 
local time of the machine A that 'corresponds' to [Tb]. It is observed that the 
25 roles of the two machines are not symmetrical and that, in general, the 
interval [Tab] is wider than [Tb]. 
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In the following, when the context does not lead to 
ambiguities, the expression 'synchronization by intervals' is quite simply 
replaced by 'synchronization'. 

By 're-synchronization' is understood more specifically a 
5 synchronization carried out after all the data has been recorded. 

The 'precision' of a synchronization by intervals is determined 
by the broadening of [Tab] with respect to [Tb]. 

The 'precision' of a synchronization process whose purpose is 
to establish a relationship between instantaneous values of Ta and Tg is the 
10 difference between the calculated value of Tab and its 'true' value. 

Calculation techniques by intervals are presented in the work 
'Applied Interval Analysis', L. Jaulin, M. Kieffer, O. Dirit and E. Walter, 
Springer-Verlag, 2001, ISBN: 1-85233-219-0. 

Figures 1A and IB illustrate how the correspondence 
15 between the times Ta and Ta, indicated by the two clocks A and B, can be 
determined. For the sake of simplicity, the discretization of the clocks has 
been neglected (6^0) and hence the lines TA(t) and TsCt), which, in reality, 
take the form of a staircase, have been approximated by straight lines. 

The first possibility, which immediately comes to mind, simply 
20 consists in reading the two clocks 'simultaneously' at a time to in order to find 
a con-espondence between TA(to) and TB(to). In fact, this 'simultaneous' 
reading is not possible in the case of real machines; moreover, this reading 
requires a finite time, which cannot be predetermined and which depends on 
the transmission times of the data over the lines or bus, together with the 

2 5 latencies of the processors and of the exploitation systems. If the demands on 

the precision of the synchronization are relatively modest, as in the case of 
multimedia applications, the imprecision in the time required to perform the 
double reading may be less than the acceptable margin of error. It is therefore 
legitimate to ignore the problem and to consider that the two readings are 

3 0 carried out substantially at the same time. 

This assumption is not justified if a higher precision is 
required. Since it is not possible to establish a one-to-one correspondence 
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Ta^Tb, an approach 'by intervals' is therefore adopted, as previously 
explained. One major advantage of this approach with respect to any 
technique that tries to determine an approximate equivalence between 
'instantaneous' dates is that it allows the uncertainty in the synchronization, 
5 which is equal to the width of the interval, to be known exactly. 

As illustrated in figure 1B, in order to determine a 
correspondence by intervals between the local times of the machines A and 
B, firstly, at a time ti, the local time of the machine A, T'a, subsequently, at a 
time to, that of the machine B, Tb, and lastly, at a time ta, again that of the 

10 machine A, T"a (it should not be forgotten that the exact values of to, ti and ta 
cannot be known) is read. It can readily be understood that it must be verified 
that no rupture in monotonicity of the local time of the machine A has 
occurred between the first and the last reading. 

One example of a method according to the invention for 

15 determining a correspondence between the times indicated by separate 
clocks and thus synchronizing data is now described in detail by means of 
figures 2, 3A - 3L, 4 and 5, together with tables Tab.A, Tab.B and Tab.C. For 
the moment, only a 'retroactive' synchronization (re-synchronization) will be 
considered, in other words where the data values are synchronized after the 

20 termination of the operation session of the system (also called 'data 
acquisition session' in the following). The case of synchronization in 'real time' 
(during the operation session of the system) will be considered later on. 

An information system, shown in figure 2, comprises three 
machines, A, B and C, such as computers, each having an internal clock Ha, 

25 Hb and He that measures a local time 7"A(t), rB(t) and 7c(t), respectively, 
together with a synchronous bus BUS, such as a FireWire® bus (IEEE 1394) 
having its own clock Hbus with a precision better than Sps, in the worst case, 
according to the IEEE1394a standard . 

The machines can be connected and disconnected from the 

30 bus BUS independently from one another; each time a reconfiguration of the 
system occurs (for example, each time a machine is connected or 
disconnected), an interrupt signal is transmitted over the bus. One of the 
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machines connected to the network is designated as 'reference machine' and 
its designation is known by the other machines; in the figure, this is machine 
B, characterized by a bold line contour. The reference machine cannot be 
considered as a 'master' machine in a centralized system: any machine in the 
5 network may be chosen to temporarily assume this function. Indeed, a new 
reference machine is chosen at random at each interrupt, which allows the 
operation of the system to be ensured even if the previous reference machine 
has been disconnected. 

Each machine is also equipped with a bus interface U, Ib and 

10 Ic, comprising a clock, HU, HIb and HIc, respectively, controlled by Hbus- The 
bus clock Hbus is generated by one of the interfaces I (called 'Cycle Master' in 
the IEEE FireWire standard). The clocks HI are feedback controlled by Hbus. 

The bus clock Hbus has a counter with limited capacity, which 
is reset to zero every 128 s, but the interface clocks comprise a software 

15 extension having a capacity that is sufficiently large so that over-runs do not 
occur over the operating duration of the system. The time measured by the 
interface clock HU, HIb and HIc ('interface time') of the machine A, B and C is 
indicated by NTA(t), NTB(t) and NTc(t), respectively. 

One machine is chosen to 'provide the network time': this 

20 means that the interface clocks of all the other machines are synchronized to 
its clock. At each inten-upt, a new machine is chosen for this task and the 
software extension of the interface clocks is reset to zero: hence, after the 
interrupt, all the clocks indicate the same time as the bus clock Hbus- 

It should be noted that the machine providing the network 

25 time is not necessarily the reference machine. If a machine is disconnected 
from the network, it is its own reference machine and its network time is the 
time on its own interface clock. 

It will be understood that the time NT of the network forms a 
first temporal reference common to all the machines connected to the bus 

30 BUS, but is non-monotonic. The time of the reference machine forms a 
second temporal reference common to all the machines connected to the bus 
BUS and, moreover, is monotonic. This second temporal reference is 
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nevertheless very imprecise because, in order to l<now it, a maclilne otiier 
tlian tlie reference one must make a request via tlie bus, whicii requires a 
time tiiat is relatively long and which cannot, in principle, be predetermined. 

Throughout its operation, each machine records, 
5 independently of the others, data values and their date stamp in the form of 
intervals, for example sensor acquisitions. 

At the same time, each machine constructs a file called 'date 
stamp file' which is formed by a series of rows. The rows are filled out at 
approximately regular intervals, for example of around one second, and they 
10 comprise five fields: 

TL1 is a first reading of the local time of the machine; 

TR is a reading of the network time; 

TL2 is a second reading of the local time of the machine; 

NMR is the designation of the reference machine; and 
15 TMR is the local time of the reference machine. 

The first three readings must be performed in order such that 
TR falls between TL1 and TL2. 

When an interrupt signal is received, the machine fills out an 
'exception row' by assigning to the fields TL1 , TR, TL2 and TMR the value 0 
2 0 and to the field NMR the conventional name 'BUS_RESET'. 

The filling out of date stamp files is a low-priority process and 
relatively undemanding, which does not interfere substantially with the 
recording and the date stamping of the events. 

Once the session is finished, the date stamp files are 
25 transferred to one and the same computer and the re-synchronization sub- 
process itself can be started. 

Here, as an example, the specific problem of re-synchronizing 
the data recorded by the machines A and B is posed. 

The tables Tab.A, Tab.B and Tab.C show a part of the date 
30 stamp file contents of the machines A, B and C, respectively, henceforth 
called FHA, FHB and FHC, respectively. Only the two rows that precede and 
follow each interrupt are shown; in addition, for the sake of clarity and saving 
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space, the field TL2 does not contain the value of the second reading of the 
local time of the machine, but the difference between the second and the first 
reading, which is clearly equivalent. The times indicated are in microseconds 
(|js), but this is not essential to the comprehension of the method. 

The interrupts subdivide each file into sections a1 - a8, pi - 
p9 and y1 - y9, separated by exception rows. It is important to observe that 
there is an interrupt at each rupture in monotonicity, hence it can be 
considered that the network time is in fact 'locally' monotonic over each 
section. 

The fields NMR and TMR allow a temporal correspondence to 
be determined between the sections of the different files. For example, it can 
be seen that the sections a1 and y2 have been recorded simultaneously, 
while the machines A and C were both connected to the bus BUS, with the 
machine A as reference machine. In this fashion, the 'history' of the system 
can be reconstructed, as illustrated in figures 3A - 3L. 

Initially, (fig. 3A), only the machine C is in operation and it fills 
out the section y1 of its date stamp file. 

Subsequently, (fig. SB), the machine A starts up and the two 
machines connect to the bus for a certain time (sections a1 and y2), A being 
the reference machine. 

Then, figure 3C, they separate again (sections a2 and y3). 

During this time (fig. 3D) the machine B goes into operation, 
but remains disconnected from the network (section p1), hence does not 
generate an interrupt. 

During the fifth period (figure 3E, sections a2, p2 and j4) the 
machines B and C are connected to the bus, and C is the reference machine. 
The machine A does not receive the interrupt signal caused by the connection 
of the machine B, hence it continues to fill out the section a2 of its date stamp 
file. 

During the sixth period, figure 3F, the three machines are 
disconnected from the bus and fill out the sections a2, p3 and y5 of the 
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respective date stamp files (the macliine A, which has already been 
disconnected, does not receive the interrupt signal caused by the 
disconnection of B and C). 

In figure 3G, the machines A and B again connect onto the 
bus, with A as reference machine, and thus pass on to the sections a3 and 
p4, whereas the machine C remains disconnected and, since it does not 
receive an interrupt signal, continues to fill out the section yS. 

Subsequently, figure 3H, the three machines are all again 
connected to the bus (sections a4, p5 and y6), with B as reference machine. 

When the latter (B) is disconnected (figure 31, sections a5, p6 
and y7), A becomes the reference machine. A short time later (figure 4J, 
sections a6, p7 and y8), B reconnects and again becomes reference machine. 

In figure 3K, the machine C disconnects (sections a7, p8 and 
y9) and A again becomes reference machine. 

Finally (figure 3L, sections a8 and p9), the machines A and B 
are also disconnected, after which the history of the system becomes 
irrelevant. 

Thanks to the second temporal reference (TMR), a kind of re- 
synchronization has thus been performed, which is however very rudimentary 
because the rows of the date stamp files are filled out at a slow rate. This 
bringing into correspondence of the sections a, p and y allows the non- 
monotonicity of the first temporal reference (NT) to be overcome and allows it 
to be used to carry out the actual re-synchronization by intervals. 

In this example, the case where only one bus is present has 
been considered. 

Coming back to the problem of re-synchronizing the data of 
the machine B with respect to the time of the machine A, the 'common 
temporal ranges' must first be determined, in other words the periods during 
which these two machines were connected to the same network via the bus 
BUS, which is that of the only network considered in this example. There are 
four of these ranges: P1 , which corresponds to the sections a3 and p4; P2, 
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which corresponds to sections a4 and p5; P3, which corresponds to the 
sections a6 and p7 and P4, which corresponds to the sections a7 and p8. 
The ranges P1 and P2, and also P3 and P4, are adjacent to one another, but 
are separated by an interrupt. 
5 Figure 4 is a graphical representation of the information 

contained in the files FHA and FHB (the discretization of the clocks Is 
neglected: the steps corresponding to the discrete increments of the counters 
are not visible). The ordinate axis carries the local times of the machines A 
and B and the time of the networks to which said machines are connected, 

10 that of the abscissae represents true time. The rows TMA and TMB represent, 
respectively, the local times of A and of B. The rows TRA (dashed) and TRB 
(circles) represent the time of the network of A and of B, respectively. The 
common ranges P1 , P2, P3 and P4 are indicated by the shaded regions. 

As was explained hereinabove, each data value recorded by 

15 the machine B is characterized by a date stamp [Tb], which is an Interval 
Included between two readings of the clock Hb surrounding a value of the 
network time, [NT]. The objective is to determine [Tab], the interval of the local 
time of the machine A corresponding to [Tb]. 

Firstly, the case where [Tb]=[Tb^] is considered, situated 

20 inside a common temporal range, for example P1, as illustrated in figure 4. 
The first operation consists in determining two date stamps [Tai^] and [Ta2^] of 
the machine A belonging to one and the same range and situated before and 
after \Tb\ In fact, these 'date stamps' are Intervals determined by the method 
described with reference to figure 1B and corresponding to the network times 

25 [NTai^] and [NTa2^], respectively. The assumption of linear drift of the clocks 
allows Uab^] to be calculated by a linear interpolation: 



[NT,' 



-[«■.,' 1 



i[T.,' 1) 



(3) 
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Since all the date stamps are intervals, the result is also an 



interval. It is clear that equation (3) could not be used if an Interrupt, and 
therefore a rupture in monotonicity of the network time, had occurred between 
[NTai^] and [NTa2^]. For this reason, it was required that [Tai^] and [Ta2^] 
belong to the same temporal range. 

The offset between the clocks of the machines A and B at the 
date [Tb] is simply given by: 



This being the difference between two intervals, the offset is 

also an interval. 

Knowing the offset between the clocks at two different dates, 
Ub^] and Ub^], the drift, assumed to be linear, can be determined: 



where [off^] and [off^] are the values of the offset between the 
clocks Ha and He at the dates [Tb^] and [Tb^], respectively. 

It is important to observe that the presence of interrupts 
between [Tb^] and [Tb^] does not constitute an obstacle for the application of 
equation (5), because the network time does not directly appear in the latter. 
On the contrary, those skilled in the art will easily understand that it is 
advantageous to maximize the separation between the dates [Tb^] and [Tb^]: 
consequently, for ]Jb\ the first date stamp of the range P1 and, for [Tb^], the 
last one of the range P4 will be taken. 

At this point, all the information required to calculate the 
correspondence between dates read on the clocks Ha and Hb, even outside of 
the common ranges P1 - P4, is available. It can indeed be shown that: 



[0ffAB'] = [TABV[TB'] 



(4) 




(5) 



lf[r, ]e [t 




then: 



[Tab]= [T^b. 



(6) 



with: 
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i^=at+|^(s=-iiB) (6.1) 

^B 

and 

'^s=n.+^^^.-^}, (6.2) 

if [Tb] <[Tb^], then: 

5 [r,, ] = [r, h[Tl ]+([r, ])-ldnft^ ];and (7) 

if Fb] >|Tb^], then: 

lT.s ]=[r, ])-[.fn/.,, ]. 8) 

[Tb]<[Tb^] is a simplified notation for <t'b and [Tb]>[Tb^] for > . 
10 Figure 5 is a graphical representation of equations 6 - 8. It 

can be observed that the width of the interval [Tab], and hence the imprecision 
in the synchronization, is more or less constant for [ r ^ ]<= |^ r ^ ,T ^ J 

and increases as it gets further away from this range (intervals [Tb'], [Tab'] and 
[Tb"], [Tab"], situated, respectively, before [Tb^] and after [Tb^]). The 
15 advantage of choosing the intervals [Tb^] and [Tb^] the furthest apart possible 
can therefore be understood. On the axes Ta and Tb of figure 5, the known 
times (i;,f;,T^,T^,l|,Tj,l|,T^,T^,i^,T^ and T^)are indicated by 
a full circle and the unknown times (T^b .T^-Ik > ^^>Z1b and T^) by a 
dashed circle. 

20 It is interesting to note that the file FHC, after having 

contributed to bringing the sections a, p and y into correspondence, is no 

longer used for the re-synchronization. 

The example of a 'hot' reconfiguration of a network discussed 

hereinabove, which is incompatible with the synchronization methods of the 
25 prior art, can now be reconsidered. First of all, it is observed that the 

reconnection of a machine Y to the network comprising a machine X poses no 
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problem of monotonicity, because the clocks of the machines are never 
synchronized to one another. Moreover, since both the machine X and the 
machine Z have been connected to the same network as the machine Y for a 
part of their history, all the data recorded by these machines can be re- 
synchronized with those of Y by a method according to the invention. 
Synchronization of the data values of the machine X with those of the 
machine Y, albeit with a lower precision, can therefore be indirectly obtained. 

In figure 6, a flow diagram is illustrated for the method of re- 
synchronization by intervals of the data of the machine B with respect to the 
machine A described hereinabove. The process can be repeated for the 
synchronization of several machines. 

The first step (E1) comprises the filling out of the date stamp 
files FHA and FHB of the machines A and B, together with those of all the 
other machines of the system (C, in particular) and, In parallel, the recording 
of the locally date stamped data values. These operations are executed up to 
the end of the data recording session. 

The second step (E2) comprises the determination, by means 
of the fields NMR and TIVIR of the files FHA, FHB and FHC, of the temporal 
relationships between the different sections of these files, and also of the 
'common temporal ranges' of the machines A and B, in other words of the 
periods during which these two machines were connected to one and the 
same network. 

Subsequently, at step E3, two rows LB1 and LB2 of the file 
FHB are chosen, each one belonging to a 'common temporal range'. These 
two rows do not need to belong to the same range: as discussed previously, it 
Is preferable that the separation between these two rows be as large as 
possible. The two readings of the local time of the machine B (TL1 , TL2) 
contained in the rows LB1 and LB2 define the intervals |TB^] and \JB\ 

At step E4, two rows (LAI and LA2) of the file FHA are 
determined that belong to the same temporal range as the row LB1 of FHB 
and are recorded before and after the latter, respectively. In the same way, 
the rows LAS and LA4 that 'surround' LB2 are determined. 



More synthetically: 

[Tai']<[Tb']<[Ta2']; 

[Tai']<[Tb']<[Ta2'] 

Advantageously, LA1 and LAS are the last rows recorded 
5 before LB1 and LB2, and LA2 and LA2 are the first rows recorded after LB1 
and LB2, respectively. 

The network times [NTai^], [NTa2^], [NTai^] and [NTa2^] are 
defined as corresponding to the field TR from the rows LA1 - LA4, 
respectively: 
10 [NTAi'] = TRrowLA1 

[NTa2^] = TR row LA2 

[NTai^] = TR row LAS 

[NTa2^] = TR row LA4 

The width of these intervals is given by the discretization 5 of 
15 the network time. 

At step E5, [Tab^], [Tab% [o^ab\ [offAs'] and [driftAB] are 
calculated by interpolation, by means of equations S - 5. 

Lastly, at step E6, the re-synchronization of all the data values 
of the machine B (or just of a part of them) is carried out by interpolation or 
20 extrapolation, by means of equations 6.1 , 6.2, 7 and 8. 

A process of re-synchronization of the data according to the 
flow diagram in figure 6 has been tested experimentally by the inventors. 

During a first experiment, two computers based on an INTEL® 
Pentium IV® processor with a clock speed of 1 .8 and 2 GHz, respectively, and 
25 a RAM memory of 512 MB, equipped with the operating system Microsoft® 
Windows 2000 Professional® and a FireWire® MindReady® interface with the 
Sednet 2® API, connected together and to a digital camera UNIBIRAIN® Fire-I® 
via a FireWire® network bus were used. The computers were to record and 
date stamp the images transmitted by the camera over the network bus; the 
30 synchronization was considered to be obtained if the date stamps of the 
images recorded by the two computers consisted of overlapping intervals. The 
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width of the re-synchronization intervals thus obtained ([Tab], again taking the 
notation of the example) was 250 - 300 jjs. 

A second experiment was carried out following the same 
protocol, but using two computers based on an INTEL® Pentium III® processor 
5 with a clock speed of 800 MHz and RAM memory of 128 and 256 MB, 
respectively, equipped with Linux RTAI® operating systems and a FireWire® 
OHCI® interface connected together and to a digital camera UNIBRAIN® Fire- 
I® via a FireWire® network bus. In this case, resynchronization intervals of 30 
- 50 MS were obtained, thanks to the fact that Linux RTAI® is executed with a 

10 real-time sub-kernel, which allows a discretization of the network time of 
around 5 ps to be attained, versus 130 ps using the Sednet 2® API in the case 
of Microsoft® Windows 2000 Professional®. 

The description presented with reference to figures 3A - 6 
and to the tables Tab.A, Tab.B and Tab.C relates to an embodiment where 

15 the date stamp files are filled out by recording two readings of the clock of 
each machine and one reading of the network time falling between the two. It 
is also possible, within the scope of the invention, to perform two readings of 
the network time and one reading of the local time of each machine, that fall 
between the two: the algorithm in figure 6 is applicable mutatis mutandis. 

20 Many other variants and improvements of the method are possible, without 
straying from the scope of the present Invention. 

For example, the determination of an offset [offAB^] (equation 
4) requires the use of a pair of readings of the network time, [NTai^] and 
[NTa2^]. In fact, any pair surrounding [Tb] can be chosen, and the result is 

25 always an interval [offAB^] containing, with certainty, the 'true' value offAB^me of 
the offset (which is impossible to know). If, starting from a plurality of such 
pairs, a plurality of intervals [offAB^]', [offAB^]", etc. is determined, it is known 
with certainty that offAB\ruee[offAB^]'n[offAB^]" n... In such a manner, a 
narrower interval, and hence a higher precision, is obtained. In the same way, 

3 0 the width of the intervals which represent the drift and the various date stamps 
Tab can be reduced. The principles of this method, referred to as 'propagation 
of constraints over the intervals' are disclosed in the article by L. Jaulin, M. 
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Kieffer, O. Dirit and E. Walter cited hereinabove. Indeed, the method of 
propagation of the constraints has allowed a reduction in the interval width of 
around 20% - 40% to be obtained for the offset ([offAB]), and around 20% for 
the drift ([driftAB])- 

5 Another improvement consists in replacing the hypothesis of 

linear drifts by a hypothesis of linearity by intervals without significant 
modifications to the algorithm. 

A further variant consists in carrying out the synchronization of 
the data during the operation session of the network: at regular intervals, the 

10 synchronization is carried out using the information available, while at the 
same time continuing to acquire data and to add rows to the date stamp files 
of the various machines. After each synchronization step, the intersection of 
the intervals thus obtained with those previously determined is performed. 

In the present document, it has always been considered that 

15 each machine assigns a date stamp to data values and records them locally. 
It will however be understood that the case in which the data is locally date 
stamped by a first machine and then transmitted over the networl< in order to 
be recorded by a second machine also forms part of the scope of the 
invention. 

20 Although, in the example considered hereinabove, the choice 

was made to synchronize the data values from two machines (A and B) with 
respect to the local time of one of them (A), the choice could also have been 
made to synchronize these data values with respect to another monotonic 
temporal reference, such as the local time of the machine C. 
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0002053250 
0003069125 



0000000012 i 
0000000035 / 

0 BUS_RESET 



A 7863621392 



0000000028 A 7908308006 

0000000029 A 7909323615 
0 BUS_RESET 0 



7937760209 
7938775811 
0 



0000000033 A 7922526444 
0000000025 A 7923542076 
0 BUS_RESET 0 



0000000029 I 

0000000028 B 8071499635 

0000000029 B 8082671426 
0000000029 B 8083687070 

0 BUS RESET 0 



7966197126 
7967212732 
0 



0000000013 A 7947916494 
0000000029 A 7948932165 
0 BUS_RESET 0 



0007804000 
0008819875 
0 



0000000030 I 
0000000026 B 8102983900 
0 BUS_RESET 0 

0000000028 A 7960103794 
0000000024 A 7961119386 



8006821287 
8007836894 
0 



0 BUS_RESET 
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Tab.B 



0002126375 0000000066 B 8022331080 

0003142125 0000000042 B 8023346654 

0010330625 0000000037 B 8030534217 

0011362125 0000000048 B 8031565573 
0 0 BUS_RESET 0 

0001864375 0000000048 C 0301503876 

0002895625 0000000025 C 0302534856 

0017335250 0000000028 C 0316972351 

0018366750 0000000038 C 0318003740 
0 0 BUS_RESET 0 



8052220716 
8053251969 
8054283215 



0001644125 
0002675500 
0003706875 



0000000042 B 8052221861 
0000000052 B 8053253220 
0000000038 B 8054284434 
0 BUS_RESET 0 



0001835750 
0002867125 



8093720757 
0 



0000000033 B 8083393725 
0000000027 B 8084424600 
0 BUS_RESET 0 

0000000035 B 8086502686 
0000000088 B 8087534237 

0000000056 B 8092690667 
0000000076 B 8093721770 
0 BUS_RESET 0 

0000000031 B 8095799819 
0000000045 B 8096831121 

0000000021 B 8101987503 
0000000037 B 8103018294 
0 BUS_RESET 0 



8114533236 
8115564503 



0000000056 B 8114534120 
0000000033 B 8115565530 

0000000078 B 8126909116 
0000000029 B 8127940507 
0 BUS_RESET 0 
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Tab.C 



0001 1 23529 000000001 6 
0002124009 0000000014 



0256150754 
0257150892 
0 



0000000015 C 0256150859 
0000000015 C 0257150996 
0 BUS RESET 0 



0273269802 
0274272003 
0 



0 BUS_RESET 



C 0298284991 
C 0299285096 
BUS RESET 0 



C 0317286987 
C 0318287084 
BUS RESET 0 



0000000014 
0000000015 
0 



C 0335288653 
C 0336288737 
BUS RESET 0 



B 8083554023 



0 BUS_RESET 



0361503412 0008986952 



A 7947723296 
A 7948733336 
BUS_RESET 0 



0370583555 
0371593569 
0 



0001851003 
0002851142 



0394604775 
0395604811 
0 



0000000015 
0000000015 



0 BUS_RESET 



